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UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 



LIBRARY ROUTINE 8 k - 212 



Exponential (D.O.I. or SADOI) 

Closed, standard entry 

21 

0, 1, 2 

Maximum error: 5 x 10" 12 (see description) 

11.3 m seconds (max) (see description) 

This routine replaces x, the contents of A before 

entry by e x . The quantity x is in the range 

-1 < x < 0. The computation is performed by eroluating 

a continued fraction expansion for e x . 

The expansion used is: 
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N. Macon, "On the Computation of Exponential and Hyperbolic Function Usinc 
Continued Fractions", Jour. ACM, PP 262 - 26?, CNstobfcr 1955. iDg 
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This routine is so witjfccn that the first 5 tens of 
7 (si* 2) are used in computing e x . It Is possible 
to nidify this, however, by changing the first digit la 
instruction 1%, [i.e. the 5 in 58r OOP] to ^ 3, or 2* 
The eaeller the mafcer used in piece of 5, the shorter 
the duration of the confutation but the greater the 
error. A simplified version of the dependency of 
speed and accuracy is given by the tables belov; 



Digit 
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for -1 < x < 0. 



TABLE I 
Max. Error 



5 x 10 
5 x 10" 
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5 x 10 
5 x 10 
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taxation 
H«3 ■ sec. 

9*7 ■ see. 

8.1 n see. 

6.5 ■ sec. 



It should be noted that the error increases with x 
approaching -1. Thus f or x in the range 
-1/2 < x < ve have the following table: 
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TABLE II 

Max. Error 

5 x 10" 12 



5 x 10 
5* 10 

5 x 10 
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LOCATION 


ORDER 


MOOES 




W *(*k) 
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K5 F 


' 


— ■ 




1 


42 16L 
L5 F 
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10 2F 
50 F 




Forai x/2 
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kO F 
7J F 




Form x /k 

—A 
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kO IF 
41 2F 
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L5 171. 
40 19L 


. 
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50 2F 
S5 F 


■ 
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10 3F 
IA 19L 
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-Compute F 
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40 20L [ 








L5 IF 
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10 5F 










66 20L 
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L5 19L 
LO 18L 




Test for end of F 
"coMpulation. 
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40 19L 
LO l8L 
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32 6L 
L9 F 




1/2 - x/2 
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Qk F* 
40 20L 




1/2 - x/2 + F/2 
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50 2F 
L5 F 
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66 20L 




xA 




Sj y I 


1/2 - x/2 + F/2 + i/2 
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00 IF 








22 F I p X 
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CRD1R 


warn 
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58? 
OOF 




Constant - ll/l6 
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Constant » l/l6 
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00 F 
OOF 
OOF 




-Temporary locations 
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